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Abstract 

Convergence results are shown for full discretizations of quasilinear parabolic partial differential 
equations on evolving surfaces. As a semidiscretization in space the evolving surface finite element 
method is considered, using a regularity result of a generalized Ritz map, optimal order error estimates 
for the spatial discretization is shown. Combining this with the stability results for Runge-Kutta 
and BDF time integrators, we obtain convergence results for the fully discrete problems. 
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1 Introduction 

In this paper we show convergence of full discretizations of quasilinear parabolic partial differential 
equations on evolving surfaces. As a spatial discretization we consider the evolving surface finite element 
method. The resulting system of ordinary differential equations is discretized, either with an algebraically 
stable Runge-Kutta method, or with an implicit or linearly implicit backward differentiation formulae. 

To our knowledge IER15j is the only work on error analysis for nonlinear problems on evolving surfaces. 
They give semidiscrete error bounds for the Cahn-Hilliard equation. The authors are not aware of fully 
discrete error estimates published in the literature. 

We show convergence results for full discretizations of quasilinear parabolic problems on evolving 
surfaces with prescribed velocity. We prove unconditional stability and higher-order convergence results 
for Runge-Kutta and BDF methods. We show convergence as a full discretization when coupled with 
the ESFEM method as a space discretization for quasilinear problems. Similarly to the linear case the 
stability analysis is relying on energy estimates and multiplier techniques. 

First, we generalize some geometric preturbation estimates to the quasilinear setting. We define a 
generalized Ritz map for quasilinear operators, and use it to show optimal order error estimates for the 
spatial discretization. During the optimal order L 2 -error bounds of the Ritz map we will use a similar 
argument as Wheeler in |Whe73j . and elliptic regularity for evolving surfaces. A further important point 
of the analysis is the required regularity of the generalized Ritz map. This will be used together with the 
assumed Lipschitz-type estimate for the nonlinearity, analogously as in (DD701ILQ951IAL151 . 

We show stability and convergence results for the case of stiffly accurate algebraically stable implicit 
Runge-Kutta methods (having the Radau IIA methods in mind), and for an implicit and linearly implicit 
fc-step backward differentiation formulae up to order five. These results are relying on the techniques 
used in TQ951IDLM12] and Al. 1 'ii 1LMV131 . By combining the results for the spatial semidiscretization 
with stability and convergence estimates we show high-order convergence bounds for the fully discrete 
approximation. 

A starting point of the finite element approximation to (elliptic) surface partial differential equations 
is the paper of Dziuk |Dzi88| . Various convergence results for space discretizations of linear parabolic 
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problems using the evolving surface finite element method (ESFEM) were shown in IDE07al IDE 13b) . a 
fully discrete scheme was analysed in [DEI2] . These results are surveyed in | DE13a] . 

The convergence analysis of full discretizations with higher-order time integrators within the ESFEM 
setting for linear problems were shown: for algebraically stable Runge-Kutta methods in IDLM12] : for 
backward differentiation formulae (BDF) in [LMV131 . The ESFEM approach and convergence results 
were later extended to wave equations on evolving surfaces, see [LM15] , 

A unified presentation of ESFEM and time discretizations for parabolic problems and wave equations 
can be found in |Manl3l . 

A great number of real-life phenomena are modeled by nonlinear parabolic problems on evolving 
surfaces. Apart from general quasilinear problems on moving surfaces, see e.g. Example 3.5 in [DE07bl . 
more specific applications are the nonlinear models: diffusion induced grain boundary motion CFP971 
IFCE011 IHan891 IDES011IES12) ; Allen-Cahn and Cahn-Hilliard equations on evolving surfaces [CENC961 
IEG961IER151IES101 IChe02] ; modeling solid tumor growth [GGG011 IES12] ; pattern formation modeled 
by reaction-diffusion equations I.B()2, 1MB14) : image processing j JYS04] : Ginzburg-Landau model for 
superconductivity I).)() 1 . 

A number of nonlinear problems, in a general setting, were collected by Dziuk and Elliott in [ DE07a , 
lDE07bl lDE13aj . also see the references therein. A great number of nonlinear problems with numerical 
experiments were presented in the literature, see e.g. the above references, in particular lDE07al lDE07bl 
IESI21IDES0T1 . 

The paper is organized in the following way: In Section [5] we formulate our problem and detail our 
assumptions. In Section [3] we recall the evolving surface finite element method, together with some of 
its important properties and estimates. We introduce the generalized Ritz map, and show optimal order 
error estimates for the residual, using the crucial IT 1,00 regularity estimate mentioned above. Section 0] 
covers the stability results and error estimates for Runge-Kutta and for implicit and linearly implicit 
BDF methods. Section [5] is devoted to the error bounds of the semidiscrete residual, which then leads to 
error estimates for the fully discretized problem. In Section [G] we briefly discuss how our results can be 
extended to semilinear problems, and to the case where the upper and lower bounds of the elliptic part 
are depending on the norm of the solution. Numerical results are presented in Section [7] to illustrate our 
theoretical results. 


2 The problem and assumptions 


Let us consider a sufficiently smooth evolving closed hypersurface T(f) C K. m+1 (m < 2), 0 < t < T, 
which moves with a given smooth velocity v. Let d m u = dtu + v ■ Vu denote the material derivative of 
the function u, where Vr is the tangential gradient given by Vru = Vit — Vzt • vv, with unit normal v. 
We are sharing the setting of [DE07al IDE13bl . 

We consider the following quasilinear problem for u = u(x,t ): 


d*u + uVrp) • v - Vr( t ) ■ (-4(M)V r (t)w) = / 

u(., 0) = u 0 


on T(f), 
on T(0), 


(1) 


where A : R. —► R is sufficiently smooth function. 

Remark 2.1. The results of the paper can be generalized to the case of a sufficiently smooth matrix 
valued diffusion coefficient A(x,t,u) : T x T(t) —> T x T(t). The proofs are similar to the ones presented 
here, except they are more technical and lengthy, therefore they are not presented here. 


The abstract setting of this quasilinear evolving surface PDE is a suitable combination of [LQ951 
Section 1] and [ AES141 Section 2.3]: Let H(t ) and V(t) be real and separable Hilbert spaces (with norms 
lUI-FZ(t), IUIv(t)> respectively) such that V(t) is densely and continuously embedded into H(t), and the 
norm of the dual space of V(t) is denoted by ||-|| v'(t)'- The dual space of H(t ) is identified with itself, 
and the duality (., .} t between V(t)' and V(t) coincides on H(t) x V(t) with the scalar product of H(t), 
for all t € [0, T], 

The problem casts the following nonlinear operator: 

(A(u)v,w)t= / A(zt)Vrz’ ■ Vrw. 

Jr(t) 
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We assume that A satisfies the following three conditions: 

The bilinear form associated to the operator A(u) :V(t)—> V(t)' is elliptic with m > 0 


{A(u)w, w) t > m\\ w||y (t) (to G V (£)), 


( 2 ) 


uniformly in u G V (t) and for all t G [0,T]. It is bounded with Ai > 0 


{A(u)v,w)t | < A^||o||y (t )||to||y(t) (v,W G V(t)), 


(3) 


uniformly in u G V(t) and for all t G [0,T]. We further assume that there is a subset S(t) C V(t) such 
that the following Lips chit z-type estimate holds: for every <5 > 0 there exists L = L{5, (S(t))o<t<T) such 
that 

||(^(toi) - A(w 2 ))u || y(t); < <5||wi -w 2 \\ V (t) +L\\wi -w 2 \\H(t), (4) 

for u G S(t), w±, w 2 G V(t), 0 < t < T. 

The above conditions were also used to prove error estimates using energy techniques in [LQ9511DD70I , 
or more recently in IAL15I . 

The weak formulation uses Sobolev spaces on surfaces: For a sufficiently smooth surface T we define 


H\T) = { v G L 2 (T) | V r ?7 £ i 2 (T) m+1 }, 


and analogously H k (T) for k G N and W k ’ p (T) for k G N ,p G [l,oo], cf. |DE07al Section 2.1]. Finally, 
Qt = U te[ o, T] r(t) x {t} denotes the space-time manifold. 

The weak problem corresponding to 0 can be formulated by choosing the setting: V(t) = H x { T(t)) 
and H{t) = L 2 (T(t)), and the operator: 

(A(u)v,w)t= / A(rt)Vrw ■ Vrw. 

Jr(t) 

The coefficient function A : K. — > K satisfies the following conditions. 

Assumption 2.1. (a) It is bounded, and Lipschitz-bounded with constant I. 

(b) The function A(s) > m > 0 for arbitrary sGl. 

Throughout the paper we use the following subspace of V(t): 


S(t.) := S(t,r) = {«G H 2 (T(t)) \ ||u|| W 2,» (r(t)) < r}. 


Then the following proposition easily follows. 

Proposition 2.1. Under Assunwtion \2.1\ and u G S{t) (0 < t < T) the above operator A satisfies the 
conditions ©, © and © (with 5 = 0). 

Proof. The first two conditions © and © are following from (a) and (b). Condition © holds, since for 
u G <S(£), wi,w 2 G i7 1 (r(£)) and any z G H 1 { r(t)), we have 



Jr(t) ^ / 

< cl llwi - w 2 || L 2 (r(t)) r ||z||iji (r(t)) , 


where the constant l is from Assumption 12.II (a). 


□ 


Definition 2.1 (Weak form). A function u G H 1 (Qt) is called a weak solution of ©, if for almost 
every t G [0, T] 



holds for every G and u(., 0) = Uq. 
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3 Spatial semidicretization: evolving surface finite elements 

As a spatial semidiscretization we use the evolving surface finite element method introduced by Dziuk 
and Elliott in |DE07aj . We shortly recall some basic notations and definitions from . 1 )K()7a . for more 
details the reader is referred to Dziuk and Elliott |Dzi881 ITTETdhl IDE13a| . 


3.1 Basic notations 

The smooth surface T(t) is approximated by a triangulated one denoted by Th[t), whose vertives are 
sitting on the surface, given as 

r hit) = U E ®- 

E(t)eT h (t) 

We always assume that the (evolving) simplices E(t) are forming an admissible triangulation Th(t), with 
h denoting the maximum diameter. Admissible triangulations were introduced in |DE07al Section 5.1]: 
every E(t) £ Th{t ) satisfies that the inner radius ah is bounded from below by ch with c > 0, and T/,(t) is 
not a global double covering of T(f). Then the discrete tangential gradient on the discrete surface Th(t) 
is given by 

v r h (t)</> := V0 - V0 • v h v h , 

understood in a piecewise sense, with i>h denoting the normal to T/j(f) (see I DEO7a) 1. 

For every t £ [0,T] we define the finite element subspace Shit) spanned by the continuous, piecewise 
linear evolving basis functions Xji satisfying t) = &ij for all i, j = 1,2 ,... ,N, therefore 

Shit) = span{\i(., t), X 2 ( ■ j t),... ,Xn(- ,<)}• 

We interpolate the surface velocity on the discrete surface using the basis functions and denote it with 
Vh- Then the discrete material derivative is given by 

dl4>h = d t 4>h + V h ■ \7(j) h {(t> h £ Shit)). 

The key transport property derived in |DE07al Proposition 5.4], is the following 

dhXk = 0 for k = 1, 2,..., N. (6) 

The spatially discrete quasilinear problem for evolving surfaces is formulated in 

Problem 3.1 (Semidiscretization in space). Find Uh £ Shit) such that 


d_ 

d t, 


\ U h 0 h + [ A{U h )V Th U h -V Th <t> h = f U h d^ h , Wh £ S h it)) 
Jt h (t) Jr h (t) Jr h (t) 


(7) 


with the initial condition Uh(. ,0) = U® £ Shi 0) being a sufficient approximation to uq. 

3.2 The ODE system 

The ODE form of the above problem can be derived by setting 

N 

Uhi-,t) =J2 a ii t )Xji-,t) 

3=1 

into o. testing with <ph = Xj an d using the transport property ©• 

Proposition 3.1 (quasilinear ODE system). The spatially semidiscrete problem Q is equivalent to the 
following nonlinear ODE system for the vector aft) = ( a.j{t )) £ M N , collecting the nodal values ofUh{-,t) : 


^^M(t)a(t)^ + A{a{t))a{t) = 0 

o(0) = ao 


( 8 ) 
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where the evolving mass matrix M(t) and a nonlinear stiffness matrix A(a(t)) are defined as 

M{t)kj = f XjXk, A(a(t)) kj = [ A(U h )Vr h Xj ' ^r h Xk, 

Jr h (t) Jr h (t) 

for a(t) defining U h = J2jLi a j(t)Xj(-,t)- 

The proof of this proposition is analogous to the corresponding one in DL.Y112 . 

3.3 Discrete Sobolev norm estimates 

Through the paper we will work with the norm and semi-norm introduced in |DLM12] . We denote these 
discrete Sobolev-type norms as 

\z{t)\ M (t) ■ = \\Zh\\L 2 (T h (t)), N(i)U(t) := \\^T h Zh\\L 2 (T h {t))i 

for arbitrary z(t) £ R N , where Zh{-,t) = ^2f =1 z j(t)Xj{ ■ >t) t further by M(t) we mean the above 
mass matrix and by A(t) we mean the linear (but time dependent) stiffness matrix: 

Mtfkj = / V rh Xj ■ Vr h Xk- 

Jr h (t) 

A very important lemma in our analysis is the following: 

Lemma 3.1 ( (DLMl2j Lemma 4.1). There are constants fx,K (independent of h) such that 

z T (M(s) - M(t))y < (e' l(s_t) - l)|^|M(t)|j/|M(t)> 
z t (A(s) - A (t))y < (e K(s_t) - 1)|^| A (t) \y\ A(t) 

for all y,z £ and s,t £ [0, T]. 

3.4 Lifting process and approximation results 

In the following we recall the so called lift operator , which was introduced in [Dzi88j and further inves¬ 
tigated in )DE07at iDEldbl . The lift operator projects a finite element function on the discrete surface 
onto a function on the smooth surface. 

Using the oriented distance function d (' IDE07al Section 2.1]), for a continuous function rjh- T hit) —> R 
its lift is define as 

VhiPit) ■■= Vh(x,t), x £ T(f), 

where for every x £ Th(t) the value p = p{x,t) £ T(f) is uniquely defined via x = p + v(p,t)d(x,t). By 
ij~ l we mean the function whose lift is rj. 

We now recall some notions using the lifting process from )Dzi88j IDE07a] and |Manl3| . We have the 
lifted finite element space 

S l h {t):={g> h =ct> l h \ ( j )h £S h {t)}. 

By Sh we denote the quotient between the continuous and discrete surface measures, dA and d Ah, defined 
as 5h&Ah = dA. Further, we recall that 

Pr := (Sij - and Pry := (% - 

are the projections onto the tangent spaces of T and T^. Further, from [DEI3b] . we recall the notation 

Qh = T-( J - dW)PrPr h Pr(7 - dU), 

Oh 

where H ('Hij = d Xj Ui) is the (extended) Weingarten map. For these quantities we recall some results 
from lDE07al Lemma 5.1], [DE13b , Lemma 5.4] and IManl31 Lemma 6.1]. 
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Lemma 3.2. Assume that Thit) and r(f) is from the above setting , then we have the estimates: 


IM||z/»(iv(t)) < ch 2 , \\vj\\L™(T h (t)) < ch, ||1 - ^llr^r,^)) < ch 2 , 
ll«9*^IU°°(r^(t)) < ch, ||Pr - Q h \\L°o(T h (t)) < ch 2 , ||Pr(^Q /l )Pr|| i =o (r , i(t)) < ch 2 , 

with constants depending on Qt, but not on t. 

Lemma 3.3. For 1 < p < oo there exists constants c\,ci > 0 independent of t and h such that the for 
all Uh £ W 1,p (Th{t)) it holds that u l h £ W 1,p (r(f)) with the estimates 

ClHw/illw^’PCIVi)) < ll M LlllV 1 .P(r(t)) < C2\\Uh\\w 1 ’P{T h (t))- 
Proof. The proofs follows easily from the relation = Pr^(J — d'H)\/ru l h , cf. IDzi88l Lemma 3]. □ 


3.5 Bilinear forms and their estimates 

Apart from the f dependence, we use the time dependent bilinear forms defined in [DEI 3b] : for arbitrary 
z,g),f £ tt 1 (r), £ £ S(t), and their discrete analogs for Zh,(fh,fh £ Sh'- 


m(z, ip) = 
a{&z,<p) = 
g(v;z,ip) = 
Kf-,v,z,g>) = 


/ Zip 

r(t) 


‘mh(Zh,4>h) — I Zfocfh 

Jr h (t) 

I A(£)Vrz-Vrip, ah((h; Z h ,(j> h ) = f A(£h)Vr h Zh ■ ^r h 4>h, 

Jr(t) Jr h (t) 

/ (V r -v)z<p, gh{V h ;Z h ,(j> h ) = [ (V r , ■ V h )Z h <j) h , 

Jr(t) Jr h (t) 

I B(&v)V r z-Vr<P, b h (f h -,V h -,Z h ,f) h )= [ B h (&; V h )V r Z h • V r ^, 
Jr(t) Jr h (t) 


where the discrete tangential gradients are understood in a piecewise sense, and with the tensors given 
as 


= 9*(A(0) + V r • M(0 - 2 A{f)V{v), 
B h ^ h \V h )ij = d m h {A{f h )) + V r , • V h A(f h ) - 2A(f h )V h {V h ), 


with 

V ( v h = ^((Vr)iVj + (Vr)jU*),, 

V h (V h )ij = ^((VrJiCVh),- + (V rh )j(V h )i), 

for i,j = 1,2,..., m + 1. For more details see iDE13bi Lemma 2.1] (and the references in the proof), or 
IDE13al Lemma 5.2]. 

We will also use the transport lemma (note that d°Zh = d t Zh + «ftVr4 for a Zh £ S l h (t)): 

Lemma 3.4. For arbitrary ff h £ S l h {t) and Zh, fh, d*Zh, d*iph £ S l h (t) we have: 

^m(z h ,iph) = m(d*Zh, <p h ) + m(z h , d*p h ) + g(v h ; z h , ip h ), 

^a(d; z h, <Ph) = a (t, l h ', d'z h , tp h ) + a(£ l h ; z h , d'p h ) + b(f, l h ;v h -z h , <p h ), 
where Vh velocity of the surface. 

Proof. This lemma can be shown analogously as |DE13bl Lemma 4.2], therefore the proof is omitted. □ 

Versions of this lemma with continuous material derivatives, or discrete bilinear forms are also true. 
The following estimates will play a crucial role in the proofs. 
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Lemma 3.5 (Geometric perturbation errors). For any f £ <S(i), and Zh,<j>h £ Sh{t) with corresponding 
lifts Zh,<Ph £ S l h (t) we have the following bounds 

| m(z h ,ip h ) - m h (Z h ,cp h ) | < c/i 2 || 2 ; ft ,|| L 2 (r(t)) ||^|| i 2 (r(t)) , 

|«.(€; Zh, <Ph) - a h (C l ]Z h , 4> h ) | < c ^ 2 IIV||i 2 (r(*)) IIVIIi 2 (r(t))» 

\9(Vh',Zh,<Ph) ~ 9hiy h ',Z h ,(j> h ) | < ch 2 \\Zh\\L 2 (T(t))\Wh\\L 2 {T{t))i 
| b{& v h',z h ,ip h ) - b h {f~ l \Vh] Zh, 4>h)\ < c/i 2 ||Vr^|| L 2( r ( t ))||Vr^||L2(r(t))- 

Proof. The first estimate was proved in |DE13bi Lemma 5.5], while the third can be found in ILM151 
Lemma 7.5]. 

The proof of the second estimate is similar to the linear case found in [DE13a , Lemma 4.7]. Again 
using the notation from [DEl3aI : 

Qh = ^-{I~ dH)PrPi h Pr{I - dU) 

Oh 

we obtain 

A{C l )Vr h Z h ■ V r h <t> h = 6 h A(C l )QhV r z h (p, ■) ■ Vm(p, •)• (9) 

Similarly as in [DE13bl Lemma 5.5], the boundedness (Proposition [2J]) and the geometric estimate 
||Pr — Qh\\L°°{Y h ) < ch 2 provides the estimate 


a{&z h ,ip h ) - a h (£ l \ Z h , <t> h ) | 

[ A(£)VrZh • Vry>/idA — f A{f~ l )S/r h Z h ■ Vr fe ^dA ft 
Jr (t) Jr h (t) 

I A(£)VrZh ■ Vr^dA — / ShA(£~ l )QhA/rZh(p, ■) • ^rPh{p,-)&A h 
Jr(t) Jr h (t) 

f A{£)(Pr-Qh)VrZh-’Vr<PhdA 
Jrit) 


< A4ch 2 ||Vr2?i||L2( r ( t ))||Vr^||L2(r(t))- 
To prove the fourth estimate we follow ILM15I : starting with the equality 

d / A~'(r‘)Vr h Z h -V Th <l>h = — 


■/ A- l (C l )Vr h Z h -Vr h ci> h = ± [ A(Z)Q l h V r z h -V r <p h 

Jr h (t) at Jr(t) 

then the transport lemma (Lemma 13.41 above) yields 

[ A(C l )dyr h Z h -Vr h <l> h + f A~ l (C l )Vr h Z h • 

Jr h (t) J r h (t) 

+ [ B h (C l ;V h )Vr h Z h -Vr k <l> h 

Jr h (t) 

= f A(OQ l h d^ r z h ■ + f A(f)Q l h V r z h ■ 

J r(t) Jr(t) 

+ / B&v h )Q l h V T z h -V T Vh+ [ d' h (A($Q l h )S7 r z h -V r <p h . 

Jr(t) J r(t) 

Therefore using that the lift of d*Zh is d'zh, © and Lemma [3© provides 

\bh{f,~ l \Vh\Zh,(j)h) - b(^;vh',Zh,^h )I 

[ dl{A{OQh)^rz h -VvV> h + [ B^v h ){Q l h -l)V T z h -Vnp h 

Jrit) Jr(t) 


'r(t) 

< ch 2 \\VrZh\\L 2 (r(t))W^rFh\\L 2 (r(t)), 

where the last estimates follow from Lemma 15721 similarly as in jLMl5l Theorem 7.5]. 


□ 
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3.6 Interpolation estimates 

By Ih : it 1 (r(t)) —> S l h (t) we denote the finite element interpolation operator, having the error estimate 
below. 

Lemma 3.6. Form < 3, there exists a constant c > 0 independent of h andt such that for u £ H 2 (T(t)^: 

ll« - + /i||Vr(u - Ihu)\\ L 2(r(t)) < ch 2 \\u\\ H 2 ^ r ^y 

Furthermore, if u £ VL 2,00 (r(t)), it also satisfies 

||Vr(ri - Ihu)\\L°°(r{t)) < c/i||«||w 2 >°°(r(t))> 
where c> 0 is also independent of h and t. 

Proof. The first inequality was shown in |Dzi88l . The dimension restriction is especially discussed in 
lDE13al Lemma 4.3]. 

The analogue of the second estimate for a reference element were shown in [SF73 i Theorem 3.1]. Then 
using standard estimates of the reference element technique we obtain the stated result, cf. IBS08) . □ 

3.7 The Ritz map for nonlinear problems on evolving surfaces 

Ritz maps for quasilinear PDEs on stationary domains were investigated by Wheeler in Whe731 . We 
generalize this idea for the case of quasilinear evolving surface PDEs. We define a generalized Ritz map 
for quasilinear elliptic operators, for the linear case see [LM15I . 

By combining the above definitions we set the following. 

Definition 3.1 (Ritz map). For a given z £ i7 1 (r(t)) and a given function £: T(t) —> K. there is a unique 
VhZ £ Shit) such that for all <j>h £ Shift), with the corresponding lift ph = <t>h, we have 

a* h iC l -,V h z,cf h ) = a*i^z,p h ), ( 10 ) 

where a* := a + m and a^ := ah + mh, to make the forms a and ah positive definite. Then VhZ £ S l h {t) 
is defined as the lift ofVhZ, i.e. VhZ = {VhZ) 1 . 

We remind here that by we mean a function (living on the discrete surface) whose lift is £. 

The Galerkin orthogonality does not hold in this case, just up to a small defect: 

Lemma 3.7 (pseudo Galerkin orthogonality). For any given £ £ S{t) there holds, that for every z £ 
H 1 {T{t)) and p h £ S l h {t) 

\a*(t;z-V h z,ph )| < c/i 2 ||T’/ l ^|| ff i ( r( t ))||^|Ui(r(t)), (H) 

where c is independent of £, h and t. 

Proof. Using the definition of the Ritz map: 

I a* if.; z - V h z, ip h ) | = | a* h {C l ;V h z, (j> h ) - a*(£; V h z, p h )\ 

< Mch 2 \\Vhz\\ H ^(T(t))\Wh\\m(T{t))i 

where we used Lemma 13.51 □ 

Error bounds for the Ritz map and for its material derivatives 

In this section we prove error estimates for the Ritz map m and also for its material derivatives, the 
analogous results for the linear case can be found in [DE13bl Section 6], !Manl3i Section 7]. The £ 
independency of the estimates requires extra care, previous results, e.g. the ones cited above, or ILM151 
Section 8], are not applicable. 




















Theorem 3.1. The error in the Ritz map satisfies the bound, for arbitrary £ £ S(t) and 0 < t < T and 
h < ho with sufficiently small ho, 

\\z ~ 'Ph,z\\L 2 (T(t)) + h\\z — 'Phz\\H 1 (r(t)) < c ^ 2 ||“||ff 2 (r(t))- 

where the constant c is independent of £, h and t (but depends on m and A 4). 

Proof, (a) We first prove the gradient estimate. 

Starting by the ellipticity of the form a and the non-negativity of the form to, then using the estimate 
CD we have: 

m\\z - V h z ||ffi ( r(t)) < a*(£; 2 - V h z, z - V h z) 

= a*(£; 2 - V h z, z - I h z) + a*(£; 2 - V h z , hz - V h z) 

< M\\z — Vhz\\H^(V(t))\\z — Ihz\\m(T(t)) 

J rch 2 \\Vhz\\m(T(t))\\IhZ — Vhz\\m(T(t)) 

< Mch\\z — VhZ\\H^(T{t))\\z\\H' 2 {T(,t)) 

+ cti 2 (2\\z - Phz\\ 2 H i^ r ^ + \\z\\ 2 H i^ r ^+ch 2 \\z\\ 2 H 2 ^ r ^, 
using the interpolation error, and for the second term we used the estimate 
\\'Phz\\m(r(t))\\hz - Vhz\\ H i(r(t)) 

< (ll'Ph.z - z|| ff i(r(t)) + ||^||/fi(r(t))) {\\Ih.z - 2 ||//i(r(t)) + ||~ - 7\~||ff 1 (r(t))) 

< 21 |z — 'Pft^||^i(r(t)) + INIfl-i(r(t)) + c ^ 2 |l-llff 2 (r(t))- 

Now using Young’s and Cauchy-Schwarz inequality, and for sufficiently small (but £ independent) h we 
have the gradient estimate 

\\ z — 'Phz\\ 2 H HY^ < —Mch 2 \\z\\ 2 H 2( T ^y 

(b) The L 2 -estimate follows from the Aubin-Nitsche trick. Let us consider the problem 
-V r • (-A(£)V r u>) + w = z - Vh.z on T(t), 
then by elliptic theory, cf. Theorem IA. 11 we have the estimate, for the solution w £ H 2 (T(t)) 

IMIff 2 (r(i)) < c\\z - Vhz\\ L 2 (r(t)), 

where c is independent of t and £. By testing the elliptic weak problem with 2 — VhZ we have 

lk-^lli 2 (r(t)) = a*(£;z-P h z,w) 

= a*(£;z ~V h z,w - I h w) + a*(£; z - V h z, I h w) 

< M\\z — Phz\\H 1 {r(t))\\w — d/tw||_f/i(r(t)) 

+ ch 2 ||'PfiZ||.f/ 1 (r(i))IC7i w; |l-tf 1 (r(i))- 

Then the estimates of the interpolation error and combination of the above results yields 

\\z ~ 'Phz\\ L 2(r(t))-\\w\\H2(r(t)) < ||Z~'P/iZ||£ 2 ( r ( f )) < Mch 2 \\z\\ H 2 ( T ( t ))\\w\\ H 2 ( T ^), 
which completes the proof of the first assertion. □ 

We will also need the following error estimates for the material derivatives of the Ritz map. 

Theorem 3.2. The error in the material derivatives of the Ritz map satisfies the bounds, for k > 1, and 
for arbitrary £ £ S(t) and 0 < t < T and h < ho with sufficiently small ho, 

k 

II K)M(z - V h z) || i2( r ( t)) + h\\V r (dl)W(z - V h z) II i 2 (r(t)) < Mckh 2 E z\\ H 2 {r(t)) . 

3 =1 

The constant Ck > 0 is independent of £ and h (but depends on a and M). 
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Proof. The proof is a modification of IManlSl Theorem 7.3]. 

For k = 1: (a) We start by taking the time derivative of the definition of the Ritz map m, use the 
transport properties ('Lemma 13.41) . and use the definition of the Ritz map once more, we arrive at 

a*(C; d*z, ip h ) = -&(£; v h ; z, ip h ) - g(v h ; z, <p h ) 

+ a ?i(£ l 'i dhPhZ, <t>h) + bh(£ l ', Vh', VhZ, 4>h) + 9h{Vh', VhZ, 4>h)- 

Then we obtain 

a *(£;d*z - d*V h z,ip h ) = - b(£;vh;z- V h z,(p h ) - g{v h ;z- V h z,g> h ) 

+ Fi{<Ph), (12) 

where 

Fi(<Ph) = (ah(r l \FhV h z,<t> h )-a%&dZV h z,<p h )) 

+ (bh(£~ l ',Vh;'PhZ,(l)h) ~ b(£-,v h ;VhZ, <p h )) 

+ ( 9h{Vh\V h z,^h ) - g{v h -,VhZ,Vh))- 

Using the geometric estimates of Lemma 13.51 F\ can be estimated as 

\FM\ < cMh 2 (\\d'Vhz\\ jji(r(t)) + ll^^llH 1 (r(t)))ll¥ , h||j? 1 (r(t))- 

Then using dffPhZ as a test function in (HD, and using the error estimates of the Ritz map, together with 
the estimates above, with h < ho independent of £, we have 

\\dh'P h z\\m(T(t)) < Mc\\d m z\\ H i(r(t)) + Mch\\z\\ H 2 (r(t))- 

Combining all the previous estimates and using Young’s inequality, Cauchy-Schwarz inequality, for suf¬ 
ficiently small (£ independent) h < ho, we obtain 

“*(£; d* z ~ dhPhZ, g>h) < Mch(\\z\\ H HT(t)) + ^||9*2:||_f/i(r(t))) (r(*))- 

Then as in the previous proof we have 

m W h z - d*P h z\\ 2 H i (m) < a*(f;d'z - d' h V h z,d' h z - d*P h z) 

= a*(£; d' h z - d*V h z, d'z - I h d*z) + a*(£; d*z - d*V h z, I h d*z - d*V h z) 

< M\\d*z - dlVhz\\m(T(t))\\dlz - Ihd*z\\m(r(t)) 

+ A1c/i^||2||ff 2 (r(t)) + || Ihd*z — d^Vhz\\H^{v(t))- 

Then the interpolation estimates, Young’s inequality, absorption using h < ho, yields the gradient esti¬ 
mate. 

(b) The L 2 -estimate again follows from the Aubin-Nitsche trick. Let us now consider the problem 
-V r • (A(€)S7rw) + w = d*z - d*V h z on T(i), 
together with the elliptic estimate (cf. Theorem lA.il) . for the solution w £ H 2 (T{t)) 

IMI-Ff 2 (r(t)) < c\\dh z ~ d' h V h z\\ LHm) , 

again, c is independent of t and £. 

Then a similar calculation as [ DE13bl . Theorem 6.2], Man 13 . Theorem 7.3] provides the L 2 -norm 
estimate. 

For k > 1 the proof is analogous. □ 
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Regularity of the Ritz map 

The following technical result will play an important role in showing optimal bounds of the semidiscrete 
residual. 


Lemma 3.8. For m < 2, there exists a constant c > 0 independent of h and t such that for a function 
u £ fL 2 ’°°(r(t)) for all t £ [0,T], the following estimate holds 

llVrT^ullioo^)) < c||tt||wa,oo(r(t)). 

Proof. Using the triangle inequality we start to estimate as 

||V r 'P/iu||i,°o( r ( t )) < || Vr(VhU - I h u) ||l°°( rp)) + ||V r (hu - u)l|L°°( r (t)) 

+ ||Vru||ioo( r ( t )). 

The last term is harmless. The second term is estimated using Lemma T3.61 For the first term, using the 
inverse estimate, error estimates for the Ritz map and for the interpolation operator we obtain 

HVrC'P/jU — Ihu)\\L°°(r(t)) < c/i -m/2 ||Vr('P/ t u - Ihu)\\ L ^{T(t)) 

< ch ~ m / 2 (|| V r (V h u - «)|| ia(r(t)) + || V r (« - 4u)|| ia(r(t)) ) 

< ch- m / 2 h\\u\\ H 2 (m) < c||u||n"2,=o(r(t))- 

□ 


Remark 3.1. A stronger result holds, assuming that u £ TU 1 , 00 (r(f)), the bound ||Vr'P/iit||L°°(r(t)) < 
c||u||w 1 .°°(r(t)) can be shown. However, the proof is technical and requires more sophisticated arguments, 
cf. fPowj /. This enables to weaken the assumption to W 1 ’ 00 in the definition of the S(t) set. We do not 
include these results here because of their length. 


4 Time discretizations: stability 


4.1 Runge—Kutta methods 

We consider an s-stage algebraically stable implicit Runge-Kutta (R K) method for the time discretiza¬ 
tion of the ODE system ©, coming from the ESFEM space discretization of the quasilinear parabolic 
evolving surface PDE. 

In the following we extend the stability result for R K methods of ■ I ) 1 . M 12‘ Lemma 7.1], to the case 
of quasilinear problems. Apart form the properties of the ESFEM the proof is based on the energy 
estimation techniques, see Lubich and Ostermann |LQ951 Theorem 1.1], Generally on Runge-Kutta 
methods we refer to 1HW96I . 

For the convenience of the reader we recall the method: for simplicity, we assume equidistant time 
steps t n := nr, with step size r. Our results can be straightforwardly extended to the case of nonuniform 
time steps. The s-stage implicit Runge-Kutta method, defined by the given Butcher tableau 


(Ci) 

( a ij ) 


(bi) 


applied to the system ©, reads as 


for i,j = 1,2, ...,s, 


S 

“I” 7” ^ ^ dijOLrij , for Z — 1, 2, . . . , 5, 

3 =1 
s 

-^n+l^n+1 — M-n®-n “1“ ^ ^ bi&ni ? 

2=1 


where the internal stages satisfy 


0 — &ni “ 1 “ 


for i = 1, 2,..., s, 
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with M n i := M(t n + Cjr) and M n+ 1 := M(t n + 1 ). Here a n i is not a derivative but a suggestive notation. 
We recall that the fully discrete solution is Ujf = 1 a n,jXj( ■ >tn)- 

For the R K method we make the following assumptions: 

Assumption 4.1. • The method has stage order q > 1 and classical order p > q + 1. 

• The coefficient matrix (a,y) is invertible. 

• The method is algebraically stable, i.e. bj > 0 for j = 1,2,. .., s and the following matrix is positive 
semi-definite: 

(bi.aij — bjOji — bibj) i j_ 1 - 


• The method is stiffly accurate, i.e. bj = a S j, and c s = 1 for j = 1, 2,..., s. 

Instead of © , let us consider the following perturbed version of the equation: 

| ^ ( M (t)a(t)) + A(a(t))a(t) = M ( t)r{t) ^ 

[ 5(0) = 5 0 . 

The substitution of the true solution 5(f) of the perturbed problem into the R-K method, yields the 
defects A ni and 6 ni , by setting e n = a n - 5(t„), E ni = a ni - a(t n + c;t) and E ni = a ni - a(t n + 
then by subtraction the following error equations hold: 


S 

M^niEni — M n e n + t ^ ( uijE n j for i — 1, 2,..., s, 

j'=i 

s 

— M n e n + t ' biE n i < 
i= 1 


where the internal stages satisfy: 


Eni T E n i — A(^<y n i')') (Xni MniXnii for i 1, 2, . . . , S, 

with r ni := r(t„ + Ciffi. 

Now we state one of the key lemmas of this paper, which provide unconditional stability for the above 
class of Runge-Kutta methods. 

Lemma 4.1. For an s-stage implicit Runge-Kutta method satisfying Assumption \4- 1\ If the equation © 
has a solution in S(t) for 0 < t < T. Then there exists a tq > 0, such that for t < tq and t n = nr < T, 
that the error e n is bounded by 


2 

M n 


n / n— 1 s n r- 

+ T J2\ ek \ 2 ^k <c(l e o]Mo + T J2J2\\Mk*rki\\l, iki + t Y1 \~ 

k— 1 ' k= 1 2=1 k= 1 


2 

M k 


+ Ct 


n— 1 

Y^ Y^ ( 


r-l, 


ki\M ki 


fc=0 i=l 



where ||w||* t = w T (A(t) +M(t)) 1 w. The constant C is independent of h, r and n (but depends on m, 
M., L, p, k and T). 

Proof. The combination of proofs of Theorem 1.1 from [L095j and of Lemma 7.1 from (DLM12] (or 
[Manl3l Lemma 3.1]) suffices, therefore it is omitted here. To be precise, the proof of this result is 
more closely related to DI..M 12 . except the estimates involving the internal stages are more similar to 
rmflbi . J □ 

Then, using the above stability results, the error bounds are following analogously as in 1 )1.M 12. 
Theorem 8.1] (or Man 13 Theorem 5.1]). 
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Theorem 4.1. Consider the quasilinear parabolic problem having a solution in S{t) for 0 < t < T. 
Couple the evolving surface finite element method as space discretization with time discretization by an s- 
stage implicit Runge-Kutta method satisfying Assumption \4-l\ Assume that the Ritz map of the solution 
has continuous discrete material derivatives up to order q + 2. Then there exists tq > 0, independent of 
h, such that for t < tq , for the error = UJf — Vh'u( ., t n ) the following estimate holds for t n = nr < T: 

n i 

IIKIU 2 (r,( t „)) + (-Ell V r,(h)^H^(r,( t3 -)))" 

3 =1 

n— 1 s i 

<Cp h , q T q+1 +C (EE \\Rh(-,tk+CiT)\\ 2 H -n rh(tk+CiT ^ + C'||l?°|| L 2( rh ( 0 )), 

fe=0 i=1 

where the constant C is independent of h, r and n (but depends on m, A4, L, /a, k and T). Furthermore 

nT 9+2 

S„ = / E»< a » ra (n«x ■A]\\L^(T h (t))At 

Jo <=1 

rT 9+1 

+ / E ( 'Ph.u)(.,t)\\ L 2(r h (t))dt■ 

Jo l=1 


The H 1 norm of Rh is defined as 


ll^+(->0ll-f/- 1 (r h (i)) sup 

o^4> h es h (t) 


{Rh(-j t), (f h )L 2 (r h (t)) 




4.2 Backward differentiation formulae 

We apply a fc-step backward difference formula (BDF) for k < 5 as a discretization to the ODE system 
©, coming from the ESFEM space discretization of the quasilinear parabolic evolving surface PDE. Both 
implicit and linearly implicit methods are discussed. 

In the following we extend the stability result for BDF methods of ILMV13 . Lemma 4.1], to the 
case quasilinear problems. Apart from the properties of the ESFEM the proof is based on Dahlquist’s 
G-stability theory (Dah78 and on the multiplier technique of Nevanlinna and Odeh 1N081 ;. 

We recall the fc-step BDF method for © with step size r > 0: 

1 k 

- y; 6jM(t n -j)a n -j + A(a n )a n = 0, (n > k), (16) 

T 3=0 

where the coefficients of the method are given by 5(£) = ?(■*■ — Cwhile the starting 

values are ao,a\,..., cxk-i- The method is known to be 0-stable for k < 6 and have order k (for more 
details, see [HW961 Chapter V.]). 

Similarly linearly implicit method modification is, using the polynomial 7 (C) = X)j=i = ( k — 

(C-l) fc 

^ k k 

- E y^ 7 ; Qri_j)q r ,. = 0, (n>k). (17) 

T 3=0 j -1 

For more details we refer to fALlS] . 

Instead of © let us consider again the perturbed problem m By substituting the true solution 
a(t) of the perturbed problem into the BDF method (fT51) . we obtain 

1 k 

''f \ djAl(t n —j^cy n —j -f- A^Oin^cXn — d n: (n 7 kf 

T 3=0 

By introducing the error e n = a n — a(t n ), multiplying by r, and by subtraction we have the error equation 
k 

E $jM„-je n ^j + rA(a n )e n + r(A(a n ) - A(a n ))a n = rd n , (n > k). 

3=0 
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In the linearly implicit case we obtain: 

k k k k 

^SjMn-jen-j + TA^'^2'y j a n -j'je n + t(a( - Af E7j“ra-j))“n = rd n , [n > k ), 

j=0 3=1 J=1 3=1 


where d n have similar properties as d n , therefore it will be also denoted by d n . 

The stability results for BDF methods are the following. 

Lemma 4.2. For a k-step implicit or linearly implicit BDF method with k < 5 there exists a tq > 0, 
such that for t < t$ and t n = nr < T, that the error e n is bounded by 

n n 

K\ m„ + tV Ma 3 < Ct V \\dj\\l >t +C max \ ei \ 2 Mi 

j=k j=k 

where ||w||2 t =w T (A(t) + M(t))~ 1 w. The constant C is independent of h,r and n (but depends on m, 
A4, L, y,, k and T). 

Proof. The proof follows the proof of Lemma 4.1 from I.MV 1.1 (using G-stability from [Dah78l and 
multiplier techniques from |N081] '). except in those terms where the nonlinearity appears. For their 
estimates we refer to Theorem 1 in [AL15] . For linearly implicit methods we follow 1AL151 Section 6]. 
Therefore these proofs are also omitted. □ 


Again, using the above stability results, the error bounds are following analogously as in 1LMV13I , 
Theorem 5.1] (or IManl31 Theorem 5.3]). 

Theorem 4.2. Consider the quasilinear parabolic problem 0, having a solution in S{t) for 0 < t < T. 
Couple the evolving surface finite element method as space discretization with time discretization by a 
k-step implicit or linearly implicit backward difference formula of order k < 5. Assume that the Ritz map 
of the solution has continuous discrete material derivatives up to order k + 1. Then there exists to > 0, 
independent of h, such that for t < tq, for the error Eff = Uf — Vh'u{-,t n ) the following estimate holds 
for t n = nr < T: 


3 =1 

n i 

< Cp h , k T k + {rJ2\\ R h{-A3)\\ 2 H-^r h (t3))) 2 + c 


3=1 


where the constant C is independent ofh.,n and r (but depends on m, M, L, p, n and T). Furthermore 




Udlf\V h u)(.A)\\L H r h (t))dt. 


5 Error bounds for the fully discrete solutions 


We follow the approach of [LMV131 Section 5] by defining the FEM residual Rh(-,t ) = Eyli r j{i)Xj{-> t) G 
S h (t) as 

[R h( j) h = ^- fVhufh + f A(Vhu)\/r(Vhu) ■ Vrfh - [(V h u)d*(/) h , (18) 

Jr h Jr h Jr h Jr h 

where (fh S Sh(t), and the Ritz map of the true solution u is given as 

N 

V h u(.,t) = 

3=1 


The above problem is equivalent to the ODE system with the vector r(t) = ( rj(t )) £ R w : 

^(M(t)a(i)) +A(a(t))a{t) = M(t)r(t), 
which is the perturbed ODE system m 


14 


















5.1 Bound of the semidiscrete residual 

We now show the optimal second order estimate of the residual Rh- 

Theorem 5.1. Let u, the solution of the parabolic problem, be in S(t) for 0 < t < T. Then there exists 
a constant C > 0 and ho > 0, such that for all h < ho and t £ [0,T], the finite element residual Rh of 
the Ritz map is bounded as 

||-Rfc||ff-i(r h (t)) < Ch 2 . 

Proof, (a) We start by applying the discrete transport property to the residual equation (1151) 

m h (Rh,<fh ) = ^m h (VhU,ct>h) + a h {V h u\ V h u, <p h ) ~ m h (V h u,dh(l>h ) 

= m. h (d*P h u, 4> h ) + a h (P h u-,V h u,(j) h ) + gh{Y h -,V h u, (f h ). 

(b) We continue by the transport property with discrete material derivatives from Lemma 13.41 but 
for the weak form, with <p := tph = (< t>hY- 

0 = ip h ) + a(u ; u, ip h ) - m(u, d*ip h ) 

at 

= m(dhU , (p h ) + a(u; u , <p h ) + g(v h \ u, ip h ) + m(u, d^iph - d*(p h ). 

(c) Subtraction of the two equations, using the definition of the Ritz map with f = u in m, i-e. 

a* h (u~ l \V h u, <f h ) = a*(u\ u , <p h ), 

and using that 

d’ph - d*iph = (' v h ~ v) - V r <Ph 

holds, we obtain 

m h {R h ,(t>h) = m. h {dlV h u,4> h ) - m{d*u,tp h ) 

+ 9h{Vh\ Vhu , 4>h) - g{vh ; u, ip h ) 

+ a* h (V h u ; V h u , (fh) - a* h (u~ l ;V h u, <j) h ) 

+ m(u, ip h ) - m h (V h u , 4> h ) 

+ m(zi, (v h - v ) • VrTh)- 

All the pairs can be easily estimated separately as c/i 2 ||<^/j||£ 2 (r(t))) by combining the estimates of Lemma 
13.51 and Theorem rrn and except the third, and the last term. 

The term containing the velocity difference (vh — v) can be estimated, using \vh — v\ + /i|Vr(uh — f)| < 
ch 2 from [|DE13bl Lemma 5.6], as ch 2 \\ Vr</3dU 2 (r(t))- 
The nonlinear terms are rewritten as: 

a* h (V h u ; V h u , <j) h ) - a* h (u~ l \P h u, <j> h ) = a* h {V h u\ V h u , <j> h ) - a*(fP h u; P h u , tp h ) 

+ a*(fP h u\ V h u, iph) - a*(u;V h u,<Ph) 

+ a*(u; V h u , <p h ) - a* h {u~ l \V h u, (j> h ) 

For the first and the third term Lemma |3.5I provides an upper bound c/z, 2 1| Vrzp*. ||z/ 2 (r(t)) (similarly like 
before). 

Finally, using Lemma 13. (31 we obtain, similarly to (p[|l . that the second term can be bounded as 
\a*(fPhU\VhU,<fh) ~ a*(u;V h u, 4> h )\ 

= [ (A{V h u) - A(u))\7 r VhU ■ Vr^ft, 

< c£\\PhU - u||L 2 (r(t))l|Vr7 : ’^'u|| i oo( r ( t ))||Vr^||L2( r ( t )) 

< cl\\VhU - zz||_L 2 (r(t)) c r || VrWft.Ilz. 2 (r(t)) 

< c£r h 2 \\Vr<Ph\\L 2 (r(t))- 

□ 
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5.2 Error estimates for the full discretizations 

We compare the lifted fully discrete numerical solution mJJ := (Ujf) 1 with the exact solution u(.,t n ) of 
the evolving surface PDE dTJ, where UJf = X^yLr Xj(-> t)> where the vectors a n are generated by a 
Runge-Kutta or a BDF method. 

Theorem 5.2 (ESFEM and R K). Consider the evolving surface finite element method as space dis¬ 
cretization of the quasilinear parabolic problem a with time discretization by an s-stage implicit Runge- 
Kutta method satisfying Assumption o Let u be a sufficiently smooth solution of the problem, which 
satisfies u(.,t ) € S(t) (T) < t < T), and assume that the initial value is approximated as 

IK - (^ , /tM)(.,0)|| L 2 (r (o)) < C 0 h 2 . 

Then there exists ho > 0 and tq > 0, such that for h < ho and r < To, the following error estimate holds 
for t n = nr < T: 


n l 

IK - w(-dn)||i2 (r(tn)) + h(rJ2 l|V r (t 3 )K “ u (->h'))ll!=(r (tj))J C(r q+1 +h 2 ). 

3=1 

The constant C is independent of h, t and n, but depends on m, At, L , p, n and T. 

Theorem 5.3 (ESFEM and BDF). Consider the evolving surface finite element method as space dis¬ 
cretization of the quasilinear parabolic problem ©, with time discretization by a k-step implicit or linearly 
implicit backward difference formula of order k < 5. Let u be a sufficiently smooth solution of the problem, 
which satisfies u(.,t) £ S(t) (0 < t < T), and assume that the starting values are satisfying 

n max IK - (P h u)(.,U) IKr(o)) < C 0 h 2 . 

U K X 

Then there exists ho > 0 and tq > 0, such that for h < ho and r < To, the following error estimate holds 
for t n = nr < T: 


n 1 

IK - u(;tn)\\L*(r(t n )) + h (f}2 ll v r(tj)K - u (-^i))llL 2 ( r( ti ))) ' < C(r k + h 2 ). 

3=1 

The constant C is independent of h, r and n, but depends on m, At, L , p, k and T. 

Proof of Theorem, \ 5. 2U 5. A The global error is decomposed into two parts: 

- u{., t n ) = (u% - ( V h u )(., t n )^j + {[V h u){., t n ) - u{., t n )j , 

and the terms are estimated by previous results. 

The first one is estimated by our results for Runge-Kutta or BDF methods: Theorem 14.11 or 14.21 
respectively, together with the residual bound Theorem 15.11 and by the Ritz error estimates Theorem 13. II 
and 13.21 

The second term is estimated by the error estimates for the Ritz map (Theorem 13.11 and 13.211 . □ 


6 Further extensions 

Semilinear problems 

The presented results, in particular Theorem 15.21 and 15.31 can be generalized to semilinear problems. 
Convergence results for BDF method were already shown for semilinear problems in |AL15j . For the 
analogous results for Runge-Kutta methods follow [L0951 Remark 1.1], Problems fitting into this frame¬ 
work can be found in the references given in the introduction. 

The inhomogeneity f(t) in the evolving surface PDE Jl]) can be replaced by f{t,u) satisfying a local 
Lipschitz condition (similar to f]J): for every 5 > 0 there exists L = L(8, if such that 

II fit, wx) - f{t , w 2 )\\ V (ty < S\\wx - w 2 \\ V (t) + L\\wi ~ w 2 \\H(t) (0 <t<T) 


16 






holds for arbitrary W\,W 2 £ V(t) with ||rui||\/(t), ||w 2 ||y(t) < r, uniformly in t. Such a condition can be 
satisfied by using the same S set as for quasilinear problems. 

To be precise: In this case the bilinear form a(.,.) is not depending on £, it is as in [DE13bj . Section[3] 
would reduce to recall results mainly from pE07al IDE13bj . The stability estimates for the Runge-Kutta 
and BDF methods are needed to be revised in a straightforward way, cf. IL095] and | lAL15] , respectively. 
The generalized Ritz map is the one appeared in 1 ..M 13. lMan!3) together with its error bounds. The 
regularity result of the Ritz map still needed from Section 13.71 

Deteriorating constants 

In view of the cited papers, especially [LQ951 Remark 1.1], Theorem 15.21 and 15.31 have an extension to the 
situation where the constants m and AT, in d2I) and ([3]) are depending on ||u||, and allowed to deteriorate 
as ||it|| tends to infinity. Using energy estimates deteriorating constants can be handled for nonlinear 
problems. Then the constant C in Theorem 15.21 and 15.31 depends also on sup tg [ 0 T ] ||u(i)||. For instance 
the incompressible Navier-Stokes equation are fitting into this framework. 


7 Numerical experiments 

We present a numerical experiment for an evolving surface quasilinear parabolic problem discretized by 
evolving surface finite elements coupled with the backward Euler method as a time integrator. The 
fully discrete methods were implemented in DUNE-FEM [ DKNOIO] , while the initial triangulations were 
generated using DistMesh |PS04] . 

The evolving surface is given by 

Y(t) = {i £ R 3 | a(t)~ 1 x\ + x 2 + £3 — 1 = 0 }, 

where a{t) = 1 + 0.25 sin(27rf), see e.g. (DE07al IDLM121 IManl3j . The problem is considered over the 
time interval [0,1]. We consider the problem with the nonlinearity A(u) = 1 — ~e~ x A The right-hand 
side / is computed as to have u(x,t ) = e~ 6t X\X 2 as the true solution of the quasilinear problem 

f d'u + vS7 r( t ) • v - Vr( t ) ■ (A(u)S7 T ( t )u) = / on T(t), 

\ m(.,0 ) = «o on T( 0 ). 

Let ( 7 fc(t))fe=i, 2 ,...,n and (Tk)k=i, 2 ,...,n be a series of triangulations and timesteps, respectively, such 
that 2 hk ~ hk -1 and 4= Tfc_ 1 , with t\ = 0.1. By efc we denote the error corresponding to the mesh 
7fc(t) and stepsize Tfc. Then the EOCs are given as 

e’z-i/'-i _ l n ( e fc/efc— 1 ) n. . 00 
EOCk — , j (k — 2, 3, ..., n). 

ln( 2 ) 

In Table [I] we report on the EOCs, for the ESFEM coupled with backward Euler method, corresponding 
to the norms 


L°°(L 2 ) : ^max N \\u^- u{.,tn)\\mr(t n )), 

N 1/2 

L 2 (H l ) : (r^||V r(tn) «- U (.,t n ))|| i2(r((n)) ) . 

n= 1 

Figure [1] shows the errors obtained by the backward Euler method coupled with ESFEM for four 
different meshes and a series of time steps. The convergence in time can be seen (note the reference line), 
while for sufficiently small r the spatial error is dominating, in agreement with the theoretical results. 
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level 

dof 

L°°(L 2 ) 

EOCs 

L 2 ^ 1 ) 

EOCs 

1 

126 

0.07121892 

- 

0.1404349 

- 

2 

516 

0.02077452 

1.78 

0.0404614 

1.80 

3 

2070 

0.00540906 

1.94 

0.0111377 

1.86 

4 

8208 

0.00136755 

1.98 

0.0033538 

1.73 

5 

32682 

0.00034289 

2.00 

0.0011904 

1.49 


Table 1: Errors and EOCs in the L°°(L 2 ) and L 2 (i7 1 ) norms 


ESFEM ESFEM 




Figure 1: ||.||M-errors of the ESFEM and the backward Euler method at time T = 1 

Figure [2] shows the errors obtained by the three step linearly implicit BDF method coupled with 
ESFEM for five different meshes and a series of time steps. Again the results are matching with the 
theoretical ones. 


ESFEM ESFEM 




Figure 2: ||.||A-errors of the ESFEM and the 3 step linearly implicit BDF method at time T = 1 
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We note that, for this example, no significant difference appeared between the fully implicit and 
linearly implicit BDF methods. 
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Appendix A: A priori estimates 

The result presented here gives regularity result, with a t independent constant, for the elliptic problems 
appeared in the proofs of the errors in the Ritz map. 

Theorem A.l (Elliptic regularity for evolving surfaces). Let T(t) be an evolving surface, fix at £ [0, T] 
and a function £: T(t) —> R. 

(i) Let f £ H -1 (r(£)) and 

L(u) := -V r • (A(£)V r u) + u. (19) 

Then there exists a weak solution u £ i7 1 (T(f)) of the problem 

L{u) = f (20) 

with the estimate 

IM|ffi(r(i)) < c||/||jj-i(r(t)), (21) 

where the constant above is independent oft. 

(ii) Let L{u) be (flT)l) . let f £ L 2 (Y(t )) and let u £ i7 1 (T(t)) be a weak solution of EOl) . Then u is 
a strong solution of dm i.e. u solves © almost everywhere and there exists a constant c > 0 
independent of t and u such that 

IMIff 2 (r(t)) < c(||u|| Z/ 2 (r ( t )) + ||/||z. 2 (r(t)))- 

Proof. For (i): The Lax-Milgram lemma shows the existence of the weak solution u. Because the coer- 
civity and boundedness constants © and © are independent of t, the constant in EXT) also not depends 
on t. For (ii): Basically we consider pullback of the operator L to T(0), rewrite it in a local chart and 
then apply the corresponding results of [GT83] . 

By assumption there exists a diffeomorphic parametrization of our evolving surface T(t), i.e. we have 
a smooth map 

4>: r(0) x [0, T] -> R m+1 

such that 

4 * 4 : r( 0 ) —> R m+1 , $ t (®) :=$(*, t) 

is an injective immersion which is a homeomorpliism onto its image with <fo(r(0)) = F(t). Because T(0) 
is compact, there exists a finite atlas 

(Wo): u n (o) cr(o)^r)" 

\ / n= 1 
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such that tp n (U n ( 0)) C R m is bounded and a finite family of compact sets (V n ( 0))^_ with V n (0) C U n { 0), 
and (Jn=i Vn( 0) = r(0). Using the properties of the diffeomorphic parametrization the new collections, 

V n (t) := $ t (y n (0)), U n {t) := $ t (J7 n (0)), := <^(0) o 

still have the same properties. Now consider the following standard formulae of Riemannian geometry 
|Fck04| : 

V7 ur J.\ Hr +\d(h° ipnit)^ 1 ) d((p n (t) 1 ) 

V r h{x,t) = 2 ^ 9n( x >t) - g- - 

i,j =1 


dxi 


where 


9ij,niX-)t) — 


d(ip„{t) *) d(<fin(t) 1 ) 


dx i dxi 

is the first fundamental form and g]%(x,t) are entries of the inverse matrix of g n := ( gij,n )? and 

V7— 1 


where X is a smooth tangent vector field with X,- = X^ • 
forward to calculate that 


dxi 


and y/tfa := 


\/det{g n ). It is straight- 


/ \ , ^ d 2 {u o ip. n [t) x ) d(uotp n it) x ) 

(—V r -.4Vr« + u) (x) = ^ a ijin {x,t) -^ b iin (x, t)-± g- - 


»d=i 

+ C n (l,t)MO( ( 3 n (f) _1 


i= 1 


for some appropriate functions £ W 1,0 ° (U n (t )), b itn ,c n £ L°°(U n (t )) where ag tn represents a 
uniform elliptic matrix. Observe that the assumptions ©, © and d4]) implies that the function above 
can be bounded independently of t. Now [ IGT831 Theorem 8.8] states that, if uoip n it) _1 is the H 1 -weak 
solution of (EU1) . then it must be a strong solution as well. 

For the estimate in (ii) observe that [GT831 Theorem 9.11] gives us for V n (t) in particular the estimate 

\\u O ip n it)~ 1 \\ H 2 (y^t)) < c(||uo^ n (t) _1 || i 2 ( C/ / i ( t )) + ||/O (*) _ 1 ||i,2(C/4(*))), (22) 

where := <p n {t)(V n (t)'\ and U’ n := <p n (t)(U n {t)) are obviously independent of t. Thus the constant 
above is independent of t. Then Theorem 3.41 in [AF031 shows that 

Nl H*(V n (t)) < C(t)||u° ^nW _1 ||ff2 ( V^ ( t) ) < C||u0^ n (t)- 1 || ff2( y j / (t)) , 

where the constant in the middle depends continuously on t, hence the last constant is independent of 
t. A similar estimate holds for the right-hand side of (l22l) . An easy calculation finishes the proof for 
(ii). ' □ 
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